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During the reporting period the following results have been accomplished: 

(1) The Incompressible Navier-Stokes solution In the vorticlty-stream 
function formulation has been applied successfully to a number 

of single airfoils at Reynolds numbers up to 2C00. Initial 
progress to higher Reynolds numbers has also been made. Results 
were presented in the reports (1), (2), and (3) below. A Pli.D. 
dissertation on this work was accepted as noted below and the 
degree was granted. 

(2) The incompressible Navier-Stokes solution in the primitive 
variable formulation (velocity-pressure) has been programmed for 
multiple airfoils with general segment arrangements and is now 
in operation and under analysis. Initial results were presented 
in the report (3) below. 

(3) The coordinate system program has been gerneralized to allow any 
segment arrangement for multiple bodies, and a number of arrange- 
ments have been investigated. Some results were presented in the 
report (3) b'low. 

(A) The potential flow solution for multiple bodies has been tested 
by comparison with the analytic solution for a circular cylinder 
pair and found to be quite good. Some results were included in 
report (3) below. 

(5) Prog amming of the compressible Navier-Stokes solution for single 
bodies using ! >rimitive variable formulation has begun. 

During the reporting period the following conference presentations 
were made (abstracts attached): 

(1) "Numerical Solution of the Navier-Stokes Equations for Arbitrary 
Two-Dimensional Airfoils," F. C. Thames, Joe F. Thompson, and 



C. W. Mnstin, Proceedings of NASA Conference on Aerodynamic 
Analyses Requiring Advanced Computers, Langley Research Center, 
Hampton, Va., March 4-6, 1975. To be published as NASA-SP 347. 

(2) "Numerical Solutions of the Unsteady Navier-Stokes Equations for 
arbitrary Bodies Using Boundary-Fitted Curvilinear Coordinates," 
J. F. Thompson, F. C. Thames, R. L. Welker, and S. P. Shanks, 
Proceedings of Arlzona/AFOSR Symposium on Unsteady Aerodynamics, 
Unlv. of Arizona, Tucson, Arlz., March 18-20, 1975. 

(3) "Use of Numerically Generated Body-Fitted Coordinate Systems for 
Solutions of the Navier-Stokes Equations," J. F. Thompson, F. C. 
Thames, C. W. Mastin, and S. P. Shanks, Proceedings of A1AA 2nd 
Computational Fluid Dynamics Conference, Hartford, Conn., June 
19-20, 1975. 

During the reporting period the following Ph.D. dissertation was 
accepted '.abstract attached): 

"Numerical Solution of the Incompressible Navier-Stokes Equations 
about Arbitrary Two-Dimensional Bodies," Frank C. Thames, Ph.D. 
Dissertation, Mississippi State University, May 1975. 


NUMERICAL SOLUT T ' N OF THE NAVIER-STOKES EQUATIONS 
FOR ARBIT TWO-DIMENSIONAL AIRFOILS 

Joe F. Thompson and Frank C. Thames 

Department of Aerophyslcs and Aerospace Engineering 
Mississippi State University, Mississippi State, MS 39762 

Research Sponsored by NASA Langley Research Center, Contract NCR 25-001-055 

Proceedings of NASA Conference on Aerodynamic Analyses Requiring Advanced 
Computers, Langley Research Center, Hampton, Va. , March 4-6, 1975 

ABSTRACT 

A method of numerical sol ;tion of the Navier-Stokes equations for the 
flow about, arbitrary airfoils, or other bodies, is presented. This meth-J 
utilizes a numerically generated curvilinear coordinate system having a 
coordinate line coincident with the body contour. Streamline, velocity 
profiles, and pressure and force coefficients for several airfoils and an 
arbitrary rock are given. Potential flow solutions are also presented. 

The procedure is also capable of treating multi-element airfoils, and 
potential flow results are presented therefor. 
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FOR ARBITRARY BODIES USING 
BOUNDARY - F1TT iD CURVILINEAR COORDINATES 

Joe F. Thompson, Frar..» C. Thames 
Ray L. Walker, Samuel P. Shanks 

Department of Aerophyslcs and Aerospace Engineering 
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Research sponsored by Langley Research Center, NASA, Grant NCR 25-001-055 
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Proceedings of Arizona/AFOSR Symposium on Unsteady Aerodynamics, Univ. of 
Arizona, Tucson, Ariz., March 18-20, 1975 

ABSTRACT 

A numerical solution of time-dependent, two-dimensional incompres- 
sible Navier-Stokes equations that can treat the unsteady laminar flow 
about bodies of arbitrary shape, such as airfoils, as naturally as simple 
bodiei has been developed. Unsteady boundaries, such as deforming bodies 
or frt e surfaces, can also be treated. This solution is based on a method 
of autt latic numerical generation of a general curvilinear coordinate system 
with coo dinate lines coincident with all boundaries of a general multi- 
connected egion containing any number of arbitrarily shaped bodies. The 
curvilinear aordinates are generated as the solution of two elliptic partial 
differential e ations with Dirichlet boundary conditions, one coordinate 
being specified t be constant on each of the boundaries, and a distri- 
bution of the other being specified along the boundaries. No restrictions 
are placed on the shape of the boundaries, which may even be time-dependent, 
and the method is not restricted to two dimensions or single bodies. 

Coordinate lines may be concentrated as desired along the boundaries. 

Spacing of the coordinate lir^s encircling the body may be controlled by 
adjusting parameters in the partial differential equations for the coordinates. 


USE OF NUMERICALLY CKNERATED BODY-FITTED COORDINATE SYSTEMS 
FOR SOLUTION OF THE NAVI ER- STORES EQUATIONS 

J. F. Thompson, F. C. Thames, C. W. Must in, anti S. P. Shanks 

Department of Aerophysics and Aerospace Engineering 
Mississippi State University, Mississippi State, MS 39762 

Research sponsored by Langley Research Center, NASA, Grant NCR 25-001-055 
and Office of Naval Research, USN, Contract N00014-74-C-0373-P0001 

Proceedings of A1AA 2nd Computational Fluid Dynamics Conference, Hartford, 

Conn., June 19-20, 1975 

ABSTRACT 

These coordinate systems remove geometric considerations from numerical 
solutions of the Navier-Stokes equations so that the flow about bodies of 
arbitrary shape can be treated. Fields with multiple bodies or with time- 
dependent boundaries may also be treated. Several types of such systems 
and methods of control of the coordinate line spacing and rlentatlon are 
discussed. Ail computations are done on a fixed rectangular field with 
square nesh without any interpolation regardless of the shape or movement 
of the phvaical boundaries and regardless of the spacing '.•* he coordinate 
lines. Potential and viscous flow solutions are presente *r airfoils 
with attached and separated flaps and for hydrofoils near a free surface. 


ABSTRACT 


Frank C. Thames, Doctor of Philosophy, 1975 

Major: Engineering, Department of Acrophysics and Aerospace Engineering 

Title of Dissertation: Numerical Solution of the Incompressible Navier- 

Stokcs Equations About Arbitrary Two-Dimensional 
Bodies 

Directed By: Dr. Joe F. Thompson 

Pages in Dissertation: 205 Words in Abstra.*: 515 

ABSTRACT 

The purpose of this research was to develop methods capable of 
producing numerical solutions of the two-dimensional, incompressible, 
time-dependent Navier-Stokes equations about arbitrary two-dimensional 
bodies. As a preliminary step a method of automatic numerical generation 
of a general curvilinear coordinate system having a constant coordinate 
line coincident with each boundary of a general multi-connected physical 
flow region was developed. These curvilinear coordinates are generated 
as the solution cf tvo elliptic partial differential equations with 
Dlrichlct boundary conditions. One natural coordinate is specified to 
be constant along each boundary of the physical region while a distri- 
bution of the other is specified along the boundary contours, "he trans- 
formation is initiated by the selection of a generating elliptic partial 
differential system of which the natural coordinates are n soli Ion in 
the physical flow plane. This generating system is transformed by inter- 
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changing the dependent and Independent variables producing a set of two 
coupled, quasl-llnear elliptic equations for the cartesian coordinates 
as functions of the natural curvilinear coordinates. If the original 
generating elliptic system obeys a maximum principle and If the range 
of the r. tural coordinate which varies along the physical contours Is 
made the same on each boundary, then *hc transformed plane will be 
rectangular. The quasl-llnear transformation equations are solved for 
the physical coordinates In the transformed plane using SOK techniques. 
The boundaries of the physical region appear simply as input boundary 
conditions for the numerical solution in the transformed plane. No 
restrictions arc placed on the shape of the boundaries in the physical 
plane which may even be time dependent. Further, the method Is not 
restricted to two dimensions. Coordinate lines may be concentrated as 
desired along the physical contours while control over the spacing of 
the radial coordinate is adjusted by varying the generating elliptic 
system. 

Once the natural coordinates are generated for a given physical 
domain, the vorticity-stream function formulation of the Navier-Stokes 
equations are transformed to the rectangular transformed plane. The 
equations are then approximated using central differences for the space 
derivatives and a backward first-order time difference for the vorticity 
time derivative. The flow solution is carried out by an implicit 
method with SOR iteration being used to converge the elliptic space 
variation at each time step. The vorticity distribution on the body 
surface is calculated at each time step utilizing a modified multi- 
dimensional false position Iteration designed to force the tangential 
velocity component on the body surface to zero. Dirichlct conditions 


corresponding to lrrotatlonal flow were used at the remote boundary. 
Successful solutions about five different bodies were generated. Those 
included a circular cylinder, cambered and flapped Karman-Tref f tz air- 
foils, a r 5ttingon 625 airfoil, and a general arbitrary body-denoted the 
cambered rock. Reynolds numbers for the solutions varied from 200 to 
2,000 at angles of attack from zero to fifteen degrees. Data defining 
stream function and vorticlty contours, boundary layer velocity profiles, 
pressure dirtributlons, and force coefficients arc presented for all 
so 1 Ions. A mildly successful solution about a NACA 0018 airfoil at a 
Reynolds number o t 10,000 was also accomplished. 


